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Abstract 



We introduce a new ^-derivable approach for the Anderson impurity model 
in a BCS superconductor. The regime of validity of this conserving theory 
extends well beyond that of the Hartree-Fock approximation. This is the first 
generalization of the [/-perturbation theory to encompass a superconductor. 



PACS Nos: 72.15.Qm, 75.20.Hr, 75.30.Mb, and 74.25.-q 



2 



Typeset using REVTpjX 



The Anderson impurity model [1] provides one of the most versatile field-theoretic de- 
scriptions for interacting correlated electrons within condensed-matter physics [2]. The 
Anderson hamiltonian in a normal metal is 

Ha = H s + H s d + Hd + Hu , (1) 

where H s = J2k,a £ kan k a describes the electron gas, H sd = Ea^^^LA + V k* d t c ka) is the 
admixture interaction, Hd = J2aE a n a represents the rf-electron level and Tiu = Un^ni de- 
notes the Coulomb-repulsion interaction. This model describes the continuous transition of 
a nonmagnetic resonant level (for f/ C T, where V = irN(0){\V\ 2 )) to a magnetic atom 
(U 3> T), and one can consider treating either V or U perturbatively. For Y/U <C 1, in 
the magnetic Schrieffer- Wolff (SW) limit [3], the Anderson hamiltonian reduces to the s — d 
hamiltonian, i.e., the Kondo model. This limit and the magnetic-nonmagnetic transition 
have been succesfully treated within the renormalization group (RNG) program [4]. How- 
ever, simpler controlled approaches would be highly desirable due to the wide applicability 
of the Anderson model and its variants to many physical systems of interest. 

After its initial introduction to describe the nonmagnetic-magnetic transition of impu- 
rities in otherwise nonmagnetic metals, and the associated many-body Kondo phenomenol- 
ogy, the Anderson model has been extensively applied and generalized to also describe 
interacting pairs of impurity atoms in metals (the Alexander- Anderson model), valence fluc- 
tuations and heavy- fermion materials (the periodic Anderson model), chemisorption (the 
Anderson-Newns model) and charging phenomena and Coulomb blockade in quantum dots 
and quantum-dot arrays. It is also of great inherent interest to consider the Anderson model 
in a superconductor and the appropriate theoretical approaches to this problem. In partic- 
ular, the RNG approach has not been generalized to this case and the Bethe-Ansatz (BA) 
method fails to be suitable since the superconducting electronic spectrum does not fulfill 
the requirement of a linear {e k oc k) dispersion relationship, necessary for the applicability 
of the /c-state enumeration within the BA scheme. 

Yosida and Yamada [5] first pointed out that in order to obtain the single-particle Green's 
function and the local density of impurity ci-states, it is useful to study the many-body 
perturbation theory with respect to U . Thus one considers the Anderson hamiltonian as 

n A = n HF -U(n^(n l )+H' u , (2) 
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where Hhf ~ ^( n t)( n l) ^ s U P to the constant energy shift, C/(n-|-)(nj_), the unperturbed 
Hartree-Fock (HF) hamiltonian and 

H'u = U5n^ ni = Ufa - (n^)^ - (n^)) (3) 

is the Coulomb- repulsion term, treated as the perturbation, see Fig. 1. Yosida and Yamada, 
in their original paper [5], used a complicated formalism utilizing Pfafnan determinants to 
derive the ?7-perturbation theory; Yamada [6] first presented the numerically evaluated d- 
electron spectral density function using the U 2 selfenergy in the low-temperature limit. This 
approach may easily be generalized to arbitrary temperatures in normal metals [7], for which 
one may introduce the ci-electron impurity selfenergy Sjv(cj) through 

G N (u) = ((G HF )- 1 -Z N (u)y 1 . (4) 

Here the HF propagator is 

G hf {uj) = -(u-E HF - F(u)) 1 , (5) 

with F(u) = Efe \V k \ 2 G° k (uj) « -iT and G° k (uj) = (u - e k y l . Above, E HF = U + (n) is the 
Hartree-Fock energy. Here and in what follows we omit the spin indices for brevity; hence 
our expressions are valid for zero field. Using the second-order free-energy functional $^ 
in Fig. 2, one easily finds that the U 2 selfenergy may be obtained by cutting the ci-electron 
propagator line in the diagram as 

J] N = 5§ N /5G N . (6) 

Consequently, the imaginary part of the second-order contribution (in U) to the impurity 
self-energy is given by 

/dcu c duo f 
~ J ~ J duj3 5{u-ui-u 2 - u 3 )F(ui, u 2 , ujz)G'' hf {ui)G'' hf {u 2 )G'' hf {u)z) . 

(7) 

Here F(ui, u 2 , to 3 ) abbreviates the following collection of thermal occupancy factors 

F(u u cu 2 , cu 3 ) = [1 - fM] [1 - f(u 2 )] [1 - f(u 3 )] + fMf(u 2 )f(u 3 ) , (8) 
with f(co) = (e^/ T + l)^ 1 denoting the Fermi distribution function. 
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Note that the [/-perturbation theory can be derived from a free-energy functional, <3>. 
Therefore, this is a "conserving approximation" [8] for the many-body system, with positive 
definite spectral functions and with sum rules fulfilled by construction. The linear term in 
U, see Fig. 2, is the HF term, describing the motion of a ci-electron with spin a in the mean 
field produced by the ci-electron with spin —a. This mean field, or the expectation value 
(n), must be computed self consistently. The quadratic term in U describes the interaction, 
at the impurity site, of the localized spin fluctuations (LSF) represented by the particle-hole 
spin-susceptibility bubbles {GnGn)- 

Yamada [6] first showed that the [7 2 -selfenergy yields a triple-peaked structure 

for the spectral density of the impurity atom. The sharp central peak for T = obtains 
the unitary limiting value at u = 0. This approach was later generalized to arbitrary finite 
temperatures [7] and it was shown that the central zero-frequency peak has a sensitive T 
dependence. 

The Anderson model in a superconductor is given by 

7~(-a,bcs = 7~(-bcs + 7~Lsd + Hd + 7~lu , (9) 
where the BCS hamiltonian is 

Kbsc = £ k*n ka - J2(^ c h c ~ki + A*c_ fci c fcT ) . (10) 

k,a k 

This model has been discussed in the HF [9, 10] and Schrieffer- Wolff [11] limits. In the HF 
approximation for a superconductor one truncates the Coulomb-interaction term as follows: 

Un^n l -> U(n ] -}n l + U{n{)n^ + £/(d{4) rf l rf t + U(did^)d\d\ . (11) 

Here the anomalous average {d^) induced at the impurity site presents another mean field, 
in addition to (n), which is to be computed self consistently. The HF approximation to 
the Anderson model in a superconductor [9] shares the same instability problem as the HF 
approximation in the normal metal [1]: a spontaneous unphysical breaking of symmetry at 
the impurity site. Our aim is to develop an approach which is free from this HF instability 
and which enables one to go beyond the HF picture. In particular, we are interested to 
investigate the qualitatively new physical features, especially in the ci-electron density of 
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states, due to increasing electron correlations as U increases beyond 71T, where the HF 
solution no longer provides quantitatively meaningful answers. 

A generalization of the [/-perturbation theory for the Anderson model in a superconduc- 
tor has thus far not been discussed in the literature. The purpose of this Letter is to suggest 
a new, conserving, self-energy [/-perturbation expansion that is valid in a superconductor. 
We note that the RNG approach has recently been generalized to treat magnetic impurities 
in superconductors, but thus far only for the s — d model [12]. 

The matrix selfenergy expansion in the Nambu space is introduced as 



GsH = ((GWMr 1 - SsH) , (12) 

where the hat denotes matrices in the particle-hole Nambu space and the ci-electron propa- 
gator in the HF approximation is given as 

. , ( u; - E HF - F u (u) Uid^-F^y 1 

G hf {uj) — — , , (13) 

\U(dl a dt) - F 21 (oj) u + E HF - F 22 {u) ) 

which, for brevity, we denote here as: 

G " M = [m cm J ' (14) 

These propagators are illustrated graphically in Fig. 3 as lines with two arrows. 

The energy-integrated Green's function (or generalized density of states) Fs(cu) in 
Eq. (13) is 

F 8 (u)=[ p V, =E^H^ fc , (15) 

\F 21 {u>) F 22 (lu) J k 

where 

V k =[ . 16 

Above, in Eq. (15), G^(u) is the unperturbed Green's function for the bulk superconductor: 

(oo — Sk A \ 1 
A* uo + e^k) 
The selfenergy in Eq. (12) is a matrix in the Nambu space 
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En (a;) E 12 (u) 
£21 M E 22 (w) 



(18) 



which may be obtained as 



Eij = 5<S>s/5(G s ) 



(19) 



where $5 now denotes the free energy in a superconductor. 

(2) 

The second-order free-energy term in U for a superconductor, $5 , is illustrated in Fig. 4 
from which we obtain the following expression, accurate to U 2 : 



where F(u>i, u 2 , 003) is, again, given by Eq. (8). Note that the propagators Q, Q, J 7 , and T in 
the above expression may be evaluated in the HF approximation, see Eq. (13), which already 
contains the pairing ineraction T~Cbcs to infinite order in the unperturbed hamiltonian 7i°. 
Therefore, the uo integrals in Eq. (20) are rather complicated: they contain delta-function 
contributions from the bound states and also continuum contributions. Trivially, one ob- 
serves that the normal-state limit, Eq. (7), is obtained consistently from Eq. (20) when 
A — > and that the diagrams in Fig. 4 reduce to those in Fig. 2 for A — > 0. 

The free-energy diagrams in Fig. 4 now comprise of two contributions linear in U, cor- 
responding to the selfconsistent occupation- number field (n), and the induced anomalous 
average (the proximity pairing at the impurity d-orbital site) (o^df). Furthermore, there 
now occur three terms quadratic in U. Two of these terms are due to the localized spin 
fluctuations (LSF), represented by the spin-susceptibility bubble (QQ), mutually interacting 
at the impurity site and also with the induced localized pairing fluctuations (LPF), shown 
as the pairing-susceptibility bubble [TT). The third contribution arises from the induced 
mutually interacting localized pairing fluctuations. 

Preliminary numerical results [13] indicate the doubling of bound states, in comparison to 
the HF theory. In particular, the new bound states tend towards uj = for increasing U. We 
shall discuss the full numerical results in detail elsewhere. Our approach can also be readily 




UJ - Ui - u 2 - uj 3 )F(uji, U 2 , UJ3) 
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extended to other situations of interest, such as an Anderson impurity in unconventional 
superconductors [14]. In this case, the self-energy expressions are formally the same but the 
order parameter must in general be interpreted as a matrix in spin space. Also the order 
parameter A(k) for unconventional superconductors possesses less rotational symmetry than 
that in the s-wave case. This will naturally lead to a /c-dependence of the <i-electron Green's 
function Gs(k,u>), the matrix selfenergy T> s (k,u) and the bound-state spectrum below the 
/c-dependent energy-gap edge. 
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FIGURES 

FIG. 1. The local <i-electron propagator in the normal state, Gn{^), is here denoted as a line 
with an arrowhead (l.h.s.). The local Coulomb-repulsion or Anderson-Hubbard interaction term 
Un^ni is represented in the diagrams with a dashed line having two vertices (r.h.s.). 

(2) 

FIG. 2. The terms $> N m the normal state up to the second order in U. The linear term in U 
is the Hartree-Fock bubble. The second-order term in U contains two particle-hole (susceptibility) 
bubbles, describing interacting localized spin fluctuations (LSF) at the impurity site. 

FIG. 3. Elements of the localized 2x2 Green's-function matrix in the Nambu space for the 
superconducting state are represented by lines with double arrows. Here Q(u), corresponding to 
the Gn(uj) in Fig. 1, is the electron propagator, while Q{oS) is the time-reversed hole-propagator 
function; is the anomalous particle-hole Green's function and T{uj) denotes its conjugate. 

(2) 

FIG. 4. Generalization of <E> to the second order in U for the pair-correlated state, <& s , ex- 
pressed in terms of the superconducting propagators, Q, Q, JF, and T, in Fig. 3. Due to the anoma- 
lous propagators in the pair-correlated medium, there now occur two first-order Hartree-Fock terms 
linear in U and three second-order terms in U, owing to the interacting localized spin fluctuations 
(LSF), correlated localized spin and pairing fluctuations (LSPF), and interacting localized pairing 
fluctuations (LPF), respectively. 
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Figure 3: 



